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ABSTRACT 


Many processes that are observable for an entire 
population can be difficult to observe for the indi- 
viduals in the population. In a conceptually similar 
problem, a production process may be influenced 
by factors that vary in a temporal sequence during 
production, but the yield is observable only at the 
completion of the process. In this note, statistical 
procedures for disaggregation, originally developed 
by R. A. Fisher, are reintroduced with suggestions 
for how they might be applied to problems of model- 
ing effects of climate and stand structure on forest 
growth. 
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How resources are allocated among individuals 
in populations is an important problem in modeling 
ecosystems. Biophysical theory, although helpful, 
is often incomplete, leaving critical parameters to 
be determined. Available radiation, nutrients, and 
moisture can be measured with relative ease on an 
area basis. How much is available to and used by 
individual trees occupying that area is either un- 
observable or very expensive to estimate. A concep- 
tually similar problem arises when total production 
can be measured only at the end of a period span- 
ning several stages of development, each of which 
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Statistical Procedures for 
Disaggregation Applicable 
to Modeling Climatic 
Effects on Forest Growth' 


may be influenced by varying conditions of environ- 
ment. The purpose of this presentation is to reintro- 
duce some statistical procedures that are applicable 
to the problem of disaggregation of effects. 
Applications of statistical procedures for disag- 
gregating effects include estimating the effects of 
daily precipitation on seasonal yield (Fisher 1925; 
Schumacher and Day 1939; Schumacher and Meyer 
1937), allocating milling costs to log size when only 
daily production totals and costs can be observed 
(Schumacher and Jones 1940), allocating growing 
space to individual trees when only stand totals can 
be determined (Chisman and Schumacher 1940; 
Curtis 1971; Lexen 1939), or estimating contribu- 
tions of different stand types to an inventory of dead 
timber when measurements are only available for 
compartment totals (Stage unpublished). All of 
these disparate examples share a common element. 
In each, the totals are available, but the individual 
contributions are either unobservable or prohibi- 
tively expensive to measure. Furthermore, in each 
case, an additive model is intrinsically reasonable. 


HISTORICAL ROOTS 


The seminal work addressing the problem of dis- 
aggregation is a paper by R. A. Fisher in 1925 titled 
“The Influence of Rainfall on Yield of Wheat at 
Rothamsted.” From that origin, the trail leads 
to a statistical training seminar in 1936 for Forest 
Service Research staff at which Fisher was a fea- 
tured lecturer. A group photo taken at that seminar 
(fig. 1) includes at least three of the authors of ap- 
plications of Fisher’s method—F. X. Schumacher 
standing to Fisher’s left, with Bessie Day, and Bert 
Lexen immediately above them. 

The procedure, although by no means new, seems 
still to have application for constructing models 
with which to evaluate effects of a globally changing 
environment. 


More than a score of forestry experts from experiment stations of the National Forest service of the United States are here attending a meeting 
at the federal building. They came here to confer with Dr. R. A. Fisher, of London, England, Fellow of the Royal Society and professor of the 
London University. The meeting will adjourn this week. Those attending are shown above. First row, left to right: Dr. G. H. Hepting, Asheville; 
Dr. |. T. Haig, Washington, D. C.; C. L. Forsling, director of the Appalachian station here; Dr. Fisher; F. X. Schumacher, Washington, D. C.; 

A. A. Hasel, Berkeley, Calif.; and T. S. Coile, Duke University, Durham. Second row, left to right: F. P. Keen, Washington, D. C.; V. L. Harper, 
New Orleans, La.; G. R. Gevorkiantz, St. Paul, Minn.; R. A. Chapman, New Orleans, La.; B. R. Lexen, Tucson, Arizona; J. H. Buell, Asheville; 
Dr. W. H. Meyer, Seattle, Washington; and Dr. A. L. MacKinney, Asheville. Third row, left to right: E. H. Frothingham, Asheville; R. A. Price, 
Ogden, Utah; K. P. Davis, Missoula, Montana; |. H. Sims, Washington, D. C.; J. G. Osborne, New Orleans, La.; Miss Bessie Day, Berkeley, 
Calif.; R. H. Blythe, St. Paul, Minn.; and C. V. Roberts, Asheville. Back row, left to right: L. E. Chaiken, Asheville; B. H. Wilford, Asheville; 

L. |. Barrett, Asheville; Dr. C. R. Hursh, Asheville; Miss Lillian Taylor, Washington, D. C.; Miss Nancy DeLachmutt, Washington, D. C.; Miss 


Dorothy Johnson, Asheville, L. H. Reineke, New Haven, Conn.; G. L. Schnur, Philadelphia, Pa.; ©. A. Connaughton, Fort Collins, Colorado; 
and C. A. Bickford, Asheville. 


Figure 1—Group of Forest Service Research staff attending Statistical Seminar at 
which Dr. R. A. Fisher was featured lecturer. Photo courtesy of Dr. Larry Davis. 
Captions supplied by Dr. Roger Chapman (taken from a newspaper account). - 


DESCRIPTION OF PROCEDURE 


The first step in the disaggregation procedure is 
to summarize the distribution of the entities con- 
tributing to the total. Several applications have 
used the first three moments of the diameter distri- 
bution (of trees in a stand, logs on a rail car, or logs 
processed each day in a mill). When the production 
units are arbitrary subdivisions of a continuous 
scale, Fisher (1925), Schumacher and Day (1939), 
Schumacher and Meyer (1937), and Zahner and 
Stage (1966) used coefficients of orthogonal polyno- 
mials describing the sequence of values (for ex- 
ample, daily rainfall during the crop year, daily 
moisture stress, or daily temperatures). For ex- 
ample, figure 2 shows a daily trend of water deficits 
as represented by a fifth-degree orthogonal polyno- 
mial. Seasonal variation of this attribute from year 
to year would be represented by variation in the 
coefficients of the polynomial fit to the data for each 
year. Application of the procedure differs between 
temporal and spatial problems. 
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Disaggregation Into a Temporal 
Sequence 


In Fisher’s original formulation, he argued that 
the effect of rainfall on wheat yield could be repre- 
sented by a multiplier of daily rainfall that changed 
slowly and continuously as the season progressed, 
and that furthermore, the multiplier could be repre- 
sented as a low order polynomial in time. For ease 
of interpretation (and computation!), use of orthogo- 
nal terms permitted one to evaluate the linear ef- 
fect, and then effects of each successively higher 
order independent of lower order effects. With the 
advent of accurate electronic computing, interest 
in orthogonal polynomials has waned, although I 
believe loss of the sequential interpretive advan- 
tages is to be regretted. 

The second step is to estimate the coefficients in 
the model. Although Fisher and his followers did 
not describe the solution of the normal equations 
as a regression problem, the calculations can in fact 
be performed with standard regression procedures. 
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Figure 2—Daily trend of water deficits for the 1958 growing season 
at the Priest River Experimental Forest, represented by fifth-degree 
polynomial. Site stores 4 inches of available water at field capacity. 


To illustrate the method of disaggregation of 
seasonal growth into daily contributions, consider 
growth as a phenomenon that starts and is com- 
pleted in four natural intervals (as shoot flushing, 
needle elongation, bud setting, and food storage, for 
example). Then it would be natural to explain the 
growth (y) as the sum of the contributions from each 
of the intervals. And in turn, the contribution in an 
interval is proportional to the causative factor, sj in 
the ith interval. 


The model is then: 
4 
a Yi caer (1) 
i=1 


and the coefficients (aj) of proportionality would be 
estimated by a least-squares method. 

If the intervals are small relative to the stages 
of development then it is reasonable to assume that 
these coefficients of proportionality change from 
interval to interval in a smooth fashion. In that 
case, they might be represented by a low-order poly- 
nomial in time. For example, we might assume a 
second-degree polynomial: 


a; = Co + C11 + coi? 
t=1,---,n (2) 


where n is the number of time intervals and the 
cj’s are coefficients in the polynomial describing 
the changing effect through time. 

Substituting this estimate of the constants in 
the original model (1): 


n 
PS S-(c +ci+ Ci”) s; (3) 


t=1 
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For perennial species and especially trees, there 
is a pronounced ante-dependence in the growth re- 
sponse. That is, the growth of the individual in the 
year previous is, per se, predictive of the growth in 
the current year. For example, if growth is poor one 
year because of adverse weather during that year, 
then growth will tend to be poor in the following 
year as well, independent of growing conditions of 
the following year. Current growth is governed not 
only by environment but also by preconditioned phys- 
iological and anatomical characteristics that may 
have been influenced by previous growing conditions. 

Thus, the model appropriate for perennial plants 
is somewhat more complex than that described by 
Fisher, Schumacher, and their colleagues. Consider 
the problem of evaluating current growth (yj) using 
only one weather factor, say moisture stress (Sj), 
where 7 is the index of year of growth. Then the 
model including effect of previous year’s growth 


(yj-1) and the weather factors in both years is: 
yi — by yi-1 = BoS; + B3S;_; + €: (5) 


where uppercase symbols are used to represent 
vector-valued variables. B and 6 are regression 
coefficients. 

In time series terminology, this model would be 
an autoregressive, moving average (ARMA) model 
but for one important difference (Monserud 1986). 
In the usual ARMA models, the S’s on the right side 
would be random, unobservable shocks to the sys- 
tem. In our case, they are observable or separately 
calculated as in Zahner and Stage (1966). Because 
of this structure of the right-hand side, the model is, 
in the terminology of econometrics, a distributed-lag 
model. 


Disaggregation to Discrete Units 


To illustrate the disaggregation to discrete units, 
consider the estimation of growing space require- 
ments (Chisman and Schumacher 1940; Lexen 
1939). The estimate of area occupied is based on 
the hypothesis that the growing space occupied by 
a single tree (and hence its share of the site re- 
sources) could be represented as a second degree 
function of tree diameter. Then the total growing 
space, which is just plot area, is the sum of growing 
space of the n individual trees, so: 


AAO Dat eee (6) 


where 


A = total land area occupied by the n trees, aid 
d = diameter at breast height of each tree. 


The b’s were estimated by least squares. Chisman 
and Schumacher (1940) further partitioned the sum 
by growth habit and species, and Curtis (1971) re- 
placed the polynomial with a power function, but 
the same procedure applies. When applying this 
technique, the choice of plots should include stands 
of widely different structure and sizes of trees, but 
with full site utilization assured by including under- 
story plants and not measuring recently disturbed 
sites. 


CRITIQUE AND OPPORTUNITIES 
FOR APPLICATION TO GLOBAL 
CHANGE MODELING 


I believe this methodology can solve some of the 
disaggregation problems in Global Change model- 
ing. There are, however, several respects in which 
modern insight into model formulation should offer 
improvement. 


Temporal Disaggregation Applications 


First, for problems of disaggregating temporal 
effects of weather within a growing season, I would 
substitute growing degree-days for Julian date. For 
example, initiation of height increment of western 
white pine depends on spring temperatures mea- 
sured by accumulation of degree days. Therefore, 
to estimate effects of moisture stress on height in- 
crement, the relation should be more accurate if 
moisture stress was measured at the same pheno- 
logical stage of development in each year rather 
than at the same Julian date. 

Next, the independent variables representing 
weather effects should be reconsidered. Zahner and 
Stage (1966) and Baier and Robertson (1968) were 
able to show that realistic modeling of soil moisture 
stress added information beyond the explanatory 
power resident in the raw weather data from which 
the soil moisture budget is calculated. Similarly, 
further transformation of soil-moisture stress into 
an index of total photosynthesis may improve the 
explanatory power of the expression even more. 

The temporal disaggregation process could be 
applied to estimate the contribution of daily photo- 
synthesis and moisture stress for the area to annual 
increments in height and diameter and to changes 
in crown length. For this step, data collected by 
Monserud and Rehfeldt at the Intermountain Re- 
search Station's Forestry Sciences Laboratory in 
Moscow would be ideal. They have stem analysis 
data describing yearly increments of height, diam- 
eter, volume, and height to crown base for three 
coniferous species. Their sample trees are in close 
proximity to the Priest River Experimental Forest 
weather station from which observations are avail- 
able since 1911, spanning most of the life of the 
sampled trees. 

When it comes to estimation, traditional least 
squares methods should be replaced by more appro- 
priate methods of cross-spectral analysis that recog- 
nize the time-series nature of the data. In such an 
analysis, the autoregressive moving-average effects 
would also be estimated. 

With this approach, it will be possible to estimate 
the contribution of daily photosynthesis to the sea- 
sonal volume increment of the sample trees. This 
information would, in turn, be very useful in resolv- 
ing the problem of carbon allocation between tops 
and roots of trees. 


Spatial Disaggregation Applications 


Leaf Area of Individual Trees—For a more 
modern application of spatial disaggregation, we 
might be interested in allocating leaf area to indi- 
vidual trees. Estimates of stand totals can be 


obtained from remotely sensed data (Running and 
others 1986) or from radiometric measurements 
beneath the canopy (Pierce and Running 1988). 
While not as appealing as direct modeling of leaf 
area of single trees (Gholz and others 1979), the 
disaggregation procedure implicitly incorporates 
within-stand interactions and correlations of errors. 


Photosynthesis for Individual Trees—Some 
models of the photosynthetic process treat the forest 
canopy as a homogeneous, three-dimensional leaf 
(for example, Running and Coughlan 1988). Disag- 
gregation of area estimates of photosynthesis and 
respiration to estimates for individual trees would 
bridge a gap between models at the landscape or 
stand level and models used to evaluate silvicul- 
tural alternatives and to translate stand growth 
into economic terms—applications that usually re- 
quire individual tree resolution. For this purpose, 
however, the procedure should be modified to take 
advantage of new information and statistical meth- 
ods. One possibility would be to replace the polyno- 
mial of tree diameter with expressions that repre- 
sent leaf area (Geron and Ruark 1988) and which 
might incorporate additional variables such as sap- 
wood basal area. In addition, Wykoff (1990) has 
shown that growth of individual trees depends on 
competition measured as a sum of growing stock in 
trees larger than the subject tree. An expression that 
represents competitive effects on growth should also 
be applicable to estimating competition for resources. 

With the modifications suggested above, one could 
disaggregate estimates of total stand photosynthesis 
into contributions from n individual trees (per unit 
of land area) using the following model: 


TP=Sch eee pA, on ian, (7) 
where 


TP = photosynthesis estimated for total stand, 
Pj = estimated allocation of daily photosyn- 
_~ thesis to ith tree 
LAj = leaf area of ith tree estimated from 
relations such as described by Gholz and 
_~. others (1979) 
LAL; = sum of leaf-area estimates over trees in 
stand larger than the subject tree 


and the 6; are estimated by econometric “errors in 
variables” methods (Johnston 1960) that are more 
efficient when the independent variables are mea- 
sured with stated errors. 

If the seasonal dynamics of leaf-area and photo- 
synthesis are not constant, one could further postu- 
late that the b’s change with time within season as 
a low-order polynomial in growing degree-days. 
Then, the methods described above could be invoked. 


One advantage of this formulation is that it obvi- 
ates the need for (or at least supplements) approxi- 
mations of light competition by Beer’s law. Further- 
more, it implicitly encompasses differential effects 
of location within the canopy on wind movement, 
temperature gradients, and competition for soil 
moisture as well as the attenuation of light. Analy- 
ses such as these would provide a cross-check on 
strictly physiological models that represent photo- 
synthesis at the level of individual trees. Through 
comparison, improvements in both approaches 
should become apparent. 
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